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Abstract 

In a Nambu-Jona-Lasinio model supplemented with an infrared cutoff in addition 
to the ultraviolet cutoff we study the issue whether diquarks are confined when the 
model is extended beyond the rainbow-ladder approximation. The gap equation, 
obtained in a truncation scheme motivated via a nontrivial quark-gluon vertex func¬ 
tion, is solved to determine the constituent quark mass if chiral symmetry is spon¬ 
taneously broken. In a second step, the Bethe-Salpeter equations for mesons and 
diquarks beyond ladder approximation are derived, taking care to preserve Gold- 
stone’s theorem in the pion channel. While the obtained masses of pseudoscalar 
and vector mesons are only moderately shifted compared to the values in ladder 
approximation, we observe that scalar diquarks disappear from the physical spec¬ 
trum and therefore are confined. For axialvector diquarks we observe indications, 
that the same mechanism may also work, but the NJL model allows no conclusive 
answer in this channel. 
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1 Introduction 


Although QCD is believed to be the theory of strong interactions hadron properties have 
not yet been calculated from this underlying fundamental theory. Thus it proves helpful to 
develop effective theories including the characterizing features of QCD as e.g. asymptotic 
freedom, confinement and chiral symmetry including its breaking pattern. A model solely 
built on chiral symmetry and leaving aside the question of confinement and the transition 
to the perturbative region is the Nambu-Jona-Lasinio (NJL) model |l|, formulated with 
quarks as elementary degrees of freedom interacting locally. Modeling the gluon sector to 
incorporate quark confinement is done in the Global Color Model (for a recent review see 
0). Both approaches, due to their construction principles, are able to describe the lowest 
lying meson states quite reasonably. Baryons then emerge as solitons of the effective 
meson theory i, 1- However, the modeled gluon interaction also gives rise to quark- 
quark correlations, leading to a picture of baryons as diquark-quark bound states || §] 
(a hybrid model combining the soliton with the bound state picture has been developed 
in ref. [0). 

Although diquarks have many appealing phenomenological aspects (see e.g. ||) and 
are useful tools to parametrize unknown and/or complicated structures, there remains 
the question whether diquarks are realized in nature. Diquarks as constituents in a 
baryon have to be necessarily in a color antitriplet state to build a colorless baryon 
when interacting with a third quark. On the other hand, an antitriplet state is forbidden 
as an asymptotic state in a presumably confining theory like QCD. When the lowest 
lying diquark states (scalar and axialvector) are evaluated in the NJL model in ladder 
approximation they are not confined. Also in the Global Color Model, which leads to 
quark confinement if certain effective quark-quark interactions (which basically have to 
provide enough interaction strength in the infrared) are used, diquarks are predicted to 
be “observable” particles 0. 

10) it has been conjectured that this fact is due to the commonly applied 


In ref. 


rainbow-ladder approximation and diquark confinement can be obtained when the quark 
Dyson-Schwinger equation as well as the diquark Bethe-Salpeter equation are considered 
beyond rainbow and ladder approximation, respectively. While these studies have been 
performed in the Munczek-Nemirovsky model JTT]] , assuming an especially simple form 
for the gluon propagator [J. it has also been claimed that the observed feature of diquark 
confinement is generic and in particular independent of the gluon propagator. Redoing 
the calculations of JiXj] with a realistic gluon propagator (e.g. with the ones reported in 

1 In all these studies it is assumed that the quark-quark interaction is given solely by the gluon propa¬ 
gator despite the fact that this is not quite correct in a non-Abelian gauge theory. E.g. ghosts invalidate 
this identification and an infrared diverging ghost propagator drastically alters this picture p2| . 
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ref. jn|) would be a very complicated task. We therefore want to test the conjecture 
of diquark confinement in a model, which assumes a quite different but also simplified 
gluon dynamics: Whereas in the Munczek-Nemirovsky model the gluon propagator is a 
delta function in momentum space leading to a constant propagator in coordinate space, 
in the NJL model the gluon propagator is assumed to be a delta function in coordinate 
space, leading to a constant propagator in momentum space. When working with a 
nonrenormalizable theory, divergent integrals have to be regularized, leaving the cutoff 
finite. We will use the proper-time method suggested in [l4[], i.e. besides the usual UV 


cutoff also an infrared cutoff is introduced to remove unphysical quark thresholds from 
correlation functions and therefore mimicing quark confinement. 

The paper is organized as follows: In the next section we solve the gap equation 
(quark Dyson-Schwinger equation) beyond rainbow approximation; in section 3 we show 
how the kernel of the meson and diquark Bethe-Salpter equation have to be consistently 
constructed to preserve Goldstone’s theorem in the pion channel. The solution of the 
Bethe-Salpeter equations for pions (0~), vector mesons (1 + ), scalar (0 + ) and axialvector 
diquarks (1 _ ) beyond ladder approximation are presented and discussed. After sum¬ 
marizing our results we will close with an outlook. In the appendices we discuss the 
elimination of thresholds in the generalized proper-time scheme and N c —counting beyond 
ladder approximation. 


2 Gap equation beyond rainbow approximation 

In Minkowski spacef] the Dyson-Schwinger equation, determining the dressed quark prop¬ 
agator is given by 

iS - 1 (p) = i> — Yj(p) 

, Mu \a \a 

E(p) = m 0 + g 2 ij k)—^S{k)T v {p,k)—. (1) 

It connects the quark propagator S(p ) with the gluon propagator D^ v {p) and the quark- 
gluon vertex function T u (p, k). The generators of color SU(3), the Gell-Mann matrices, 
are denoted by A a /2. In the NJL model the gluon propagator is assumed to be 

9 2 D tl v(p) = Gg^„, ( 2 ) 

leading to an effective, nonrenormalizable local quark-quark interaction governed by the 
coupling strength G. When furthermore the quark-gluon vertex function is substituted 

2 All calculations are done in Minkowski space; a Wick rotation is performed only to evaluate the final 
integrals. We furthermore assume exact flavor SU(2) symmetry and suppress all explicit flavor indices. 
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Figure 1: The quark Dyson-Schwinger equation including the nontrivial quark-gluon vertex 
function is shown. Note, that due to the ansatz of the gluon propagator (eq. ffty) the gluon 
lines are momentum independent and merely represent the color flow through the diagrams. 


by its perturbative expression (q^) one ends up with the quark Dyson-Schwinger equation 
of the NJL model in mean field or rainbow approximation. As suggested in ref. |T0|, a 


minimal way to go beyond rainbow approximation is to include the vertex correction to 
first order in G: 


X a 

T 


/* \U \LL \U 

r„(*, p) = — + Gzj j^- 4 D pu (p - l)r Y S(l + k - p) lp -S(l)Y Y (3) 

( 4 ) 


A° 


7m + 


The part of the vertex function T L ({k — p) 2 ), describing the dressing of the perturbative 
quark-gluon vertex is given by 

r d A l 

7 M r L ((/c - p) 2 ) = i J uS(l + k- p)^S(l)Y- 


( 5 ) 


leading to the quark Dyson-Schwinger equation 

S(p) = m 0 + ^Gi J 

—? r d A k 

+ J ^j,S(k)rr L ((k - p ) 2 ). ( 6 ) 

Throughout the paper we will denote factors, arising from the color structure by “sans 
serif’ letters. In fig. ([!]) the quark Dyson-Schwinger equation beyond rainbow approxi¬ 
mation is shown. Note, the gluon lines appearing in the Feynman diagrams are given by 
eq. (^) and are therefore momentum independent. 

Here and in the following we regularize loop integrals in a generalized proper-time 
scheme, 

W = (PijI I exp( - sJV) (PT y.If dss ^ (7) 

where N n is the denominator of a Wick rotated loop integral after introducing the ap¬ 
propriate Feynman parameter integral and performing a corresponding momentum shift. 
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Rainbow, O(G) 

Beyond Rainbow, 0(G 2 ) 


mo = 0 

m 0 = 0.017 

mo = 0 

m 0 = 0.017 

p = 0.1 

0.3873 

0.400 

0.4194 

0.4295 

p = 0.2 

0.3873 

0.400 

0.4194 

0.4297 


Table 1: The solution of the gap equation in rainbow (0(G)) and beyond rainbow approx¬ 
imation (0(G 2 )) is shown. We used parameters, which are fixed in the standard NJL 
model: The chosen set (A = 0.630 GeV, G = 184.18 GeV -2 , m 0 = 0.017 GeVj lead to 
M = 0.4 GeV and to the pionic observables m T = 0.140 GeV and f n = 0.093 GeV, when 
the dominant amplitude of the pion (see eq. 0 is considered) W- In the chiral limit 
the same parameters are used, but m 0 = 0. All masses and scales are given in units of 
GeV. 


As has been shown in ref. |T4| any nonvanishing infrared cutoff p (with p <C A, where 
A is the ultraviolet cutoff) removes unphysical quark thresholds from correlation func¬ 
tions and therefore simulates quark confinement in a crude way, for details we refer to 
Appendix A. Note, however, although this prescription allows to investigate correlation 
functions at external momenta P 2 > 4 M 2 , i.e. above the “pseudo”threshold, without 
opening the unphysical decay channel, this quark confinement mechanism is certainly too 
naive to explain physics far above the pseudothreshold. Due to the gauge invariance of 
the proper-time regularization, we especially find that the vertex correction is purely lon¬ 
gitudinal (see eq. (f|)). We furthermore observe that the vertex correction can be very 
well approximated by a simple dipole 


T l (Q 2 ) ~ 


a 

W+W' 


( 8 ) 


with coefficients a(M, A, p) and b(M,A,p) obtained by a y 2 —fit to the regularized result 
of eq. (j5|) which obviously depend on A, p and M. Using the effective parametrization of 
the quark-gluon vertex correction (|8|) in eq. (p[) one finally ends up with the gap equation 
determining the constituent quark mass Q: 

M = m 0 + ■4M 2 G(r(-l,M 2 /A 2 )-T(-l,M 2 / / r 2 )) 

+ ^7 ' If ' AMa ° 2 ( J Q dx ( l ~ ^)^( ex P(- y 7 A2 ) - exp(-V 2 /^ 2 ))) , 
V 2 = (1-z) -b + xM 2 . (9) 


While the first line of eq.([|) is the familiar mean field result, the second line arises from 
the vertex correction. Although the term induced by the vertex correction leads to a 
3 r(n, x) denotes the incomplete Gamma function arising in the proper-time regularization scheme. 
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Figure 2: In this figure we display how the kernel of the meson Bethe-Salpeter equation is 
obtained from the quark self-energy diagrams by cutting all possible internal quark lines. 
See also fig. (1) of ref. ffUf. 


momentum dependent quark mass we assume that the momentum dependence at low 
energies is weak and it is therefore legitimate to work with £(p 2 ) = M = const, in eq. 
(||) . The gap equation ([|) can be solved by taking the dependence of the dipole coefficients 
a(M, A, p) and b(M, A, p) on M properly into account. Our results are displayed in table 
([D. It is seen, that the vertex correction is repulsive, leading to slightly larger constituent 
quark masses. When doing the calculation with different values of the infrared cutoff /i, 
we find that the results are quite insensitive on p, as long as there is enough phase space 
available i.e., the condition p ■C A is fulfilled. 


3 Bethe-Salpeter equation for mesons and diquarks 

The Bethe-Salpeter equation describing mesons as quark-antiquark bound states is written 
as 

r d A k 

p) = J K M (k,p; P)(S(k + |P)<h M (P, k)S{k - \P )). (10) 

In ladder approximation the kernel KM(k,p] P) reduces itself to a momentum indepen¬ 
dent contact interaction, reflecting the local quark-quark interaction of the NJL model. 
It is well known, and can be proved quite easily, that the ladder approximation preserves 
Goldstone’s theorem: When the constituent quark mass obtained from the gap equation 
(eq. (^) without the vertex correction) in the chiral limit is used in eq. m , one imme¬ 
diately folds massless bound states with quantum numbers of pions (Goldstone bosons), 
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reflecting spontaneously broken chiral symmetry. Following ref. 110] we construct the 


kernel of the Bethe-Salpeter equation for mesons in a way which preserves the Goldstone 
theorem in every order of the truncation scheme: The Bethe-Salpeter kernel follows via 
the replacement |10 


—> 7 p S{k + \P)$ M (k, P)S(k - \P) 7 ^ 


( 11 ) 


from the quark self-energy E(p). As it is displayed in fig. (H), this corresponds to cut 
all internal quark lines in the quark self-energy diagrams, fig. (0), thereby generating the 
Bethe-Salpeter kernel. Since the diagram in 0(G 2 ) contains three internal quark lines, 
the Bethe-Salpeter kernel beyond ladder approximation consists of three parts. This 
prescription leads to the meson Bethe-Salpeter equation 

4 r <i 4 k —? r d A k r rl 4 l 

*m(P,p) = - z G m i\ ^ 7 uS(k + 1 P)*m(P, k)S(k - \P)Y + ~^~G 2 M i J — J — x 

[ 7 v S(k + \ P)$ M (P, k)S(k - \P) lp S(l + k- p)YS(l)Y 
+ 7 u S{k) lp S{l + k-p+ 1 P)$m{P, l + k- p)S(l + k-p- \P)YS(l)Y 
+ 7 v S{k) lp S{l + k- p)YS(l + I P)$m(P, 1)S{1 - \P)Y} . 

( 12 ) 


By Gm we collectively denote the coupling constants Gm = {G ni G PtU } in the different 
meson channels. While chiral symmetry forces the choice G n = G (see table (HD), the cou¬ 
pling constant in the vector meson channel can be different from that value. Nevertheless, 
we use Gm = G^ — G PjUJ — G. 

Note that the three terms contributing to the Bethe-Salpeter equation in 0(G 2 M ) 
appear with an equal color factor of —(N 2 — 1)/(4A(?) = —2/9, see Appendix B. However, 
after taking the relevant Dirac traces, the second term of order G 2 M gets an additional 
minus sign, making this term repulsive compared to the first and third term which are 
attractive. 

In order to be consistent with the approximation employed in the gap equation, i.e. 
working with a momentum independent constituent quark mass, we have to neglect the 
relative momentum of the Bethe-Salpeter amplitudes, thereby using <E>(P,p) = <F(P). For 
a pion the appropriate Dirac structure of the amplitude is given by 


MP) = 75^1 (i 3 ) + 75 # ’Mi 3 )- (I 3 ) 

We therefore consider not only the dominant Dirac structure (~ 75 ), but also the sub¬ 
dominant amplitude (~ ^ 75 ), which is necessary to respect chiral symmetry and which 
induces significant contributions to pionic observables (for a recent investigation see ref. 


7 













|T6| ). Physically, this is due to a possible tt — A-mixing. In case of vector mesons (p, u), 
the dominant transversal Dirac structure of the amplitude reads 


= (7, 

Equivalently one may use 


p^(p )=o. 


i„jp) = ejyhnn 


4p* = o. 


(14) 


(15) 


involving the polarization vector (A = 0, ±1) being orthogonal to the total momentum 
of the vector meson. A possible subdominant amplitude would correspond to a mixing 
between vector and tensor mesons, which are known to be absent in the NJL model. 
Therefore we have to restrict ourselves to the dominant Dirac structure but choose it 
in a way which incorporates transversality from the very beginning. After inserting the 
ansatze for the Bethe-Salpeter amplitudes (eq. ( |T3| ) and eq. ( P~I| ) or (|T^)) into eq. (|T2|), 
and taking appropriate Dirac traces, the Bethe-Salpeter equation can be reduced to a 
system of two coupled algebraic equations for the pions and to one algebraic equation 
for the vector mesons, with the generic structure of the inverse propagators of composite 
mesons 


[A(P 2 ) — l]p2 =M 2 4/ m (P) — 0, 


(16) 


where M m = M n or M p ^. The meson masses are given either by the zero of det[A(P 2 ) — 1] 


m 


(in the pion channel, where A(P 2 ) is a 2 x 2 matrix) or simply by the zero of eq. ([T 
the p,u channel, where it is just one equation). 

In Appendix A we discuss the pion Bethe-Salpeter equation in ladder approximation 
in order to demonstrate how the confinement mechanism via the infrared cutoff p works. 

The two-loop integrals appearing in eq. (|12D are calculated as follows: First we perform 
one integration (e.g. over the momentum /) and find that the result (“longitudinal” and 
“transversal” part) can be approximated by dipols, in analogy to eq. (H); the second 
integration can then be performed using the dipols as effective interaction between the 
quark and antiquark. The final expressions for the different A(P 2 ) are rather lenghty and 


will be explicitly given elsewhere [17 


In order to calculate scalar and axialvector diquarks, the partners of the pseudoscalar 
and vector mesons, respectively, the Bethe-Salpeter equation has to be modified in the 
following way. The model independent expression of a diquark Bethe-Salpeter equation 
can be written as 

*d(P,p) = j J^Vi K D( k ,p', P)(S(k + k P)t D (P, k)S T (k - i P)), 


(17) 









where “T” denotes the transposed of the Dirac matrix S(p). As has been shown in 
ref. 0 , each contribution to the diquark kernel Kjj(k,p] P) can be obtained with the 
replacement 


s (ph»Y ^ ( 18 ) 

from the meson kernel P) (eq. (| 12 ])) which leads to the diquark Bethe-Salpeter 

equation 

$ D (P,p) = -Y G d*J + \ P )®d{P, k)S T (-k + \P)Y T + G 2 D ijJ 

[ l - lv s{k + ±P)<M^ k)S T {-k + \P)rfS T (-l - k + p)Y T S T (-l)Y T 
+ ^iuS(k) lp s(i + k - p + fP)$ D (P, i + k - p)s T (—i -k + p + \P)Y T S T {-1)Y T 
+ \lvS{k) lp S(l + k - p)YS{l + \P)*d{P, l)S T (-l + |P)7 pT ] , (19) 


with Gn = {G*o+) G\G\. To match our result in the scalar diquark channel (when only the 
dominant Dirac structure is considered) with the ones obtained in a ladder NJL calculation 


T5| we choose Gq+ = 2/3 G. To simplify comparison, in the axialvector channel we use 


Gi~ = 2 G, thereby obtaining in ladder approximation the same masses for axialvector 
diquarks and vector mesons. 

When comparing the diquark Bethe-Salpeter equation with the one for mesons, the 
most striking difference are the different color factors in eq.(|l9|) arising from the projection 
onto color antitriplet states. While in ladder approximation the Bethe-Salpeter equations 
only differ by a factor 2 , leading to a weaker binding of diquarks when compared to 
mesons, beyond ladder approximation these color factors are more significant. The color 
prefactors of the first and third term in 0{G 2 D ) are given by (N c + l)/(2A r c ) 2 = 1/9, while 
for the second term one finds (N c + 1)(1 + N c — N^)/{2N C ) 2 = —5/9, see Appendix B. 
Combined with factors arising from the Dirac structure of these terms, the part ~ —5/9 
is repulsive, giving the possibility to avoid binding, thereby confining the diquarks. It 
is amusing to note, that for N c = 2, the color factors for mesons and diquarks (which 
are then colorless “baryons”) are identical up to a minus sign, which is compensated by 
the different Dirac structure. Unconhned mesons then automatically lead to unconfined 
diquarks, as expected. In the limit N c —> oc only the repulsive second term of eq. 
(P~9|) , which is of 0(N C ), survives (the color prefactor in ladder approximation is given by 
— (N c + 1)/2IV C = —2/3 , i.e. of 0(1)), which is consistent with Witten’s conjecture fT8[| , 
stating that mesons are the relevant degrees of freedom for N c oo. However, in QCD 
we have neither N c = 2 nor N c = oo, but N c = 3. 


9 














Solution of Pion/0 + -Diquark BSE 


in chiral limit m 0 =0 



P 2 /GeV 2 


Figure 3: The inverse of the propagators for the pion and the scalar diquark are shown 
in ladder ( 0(Gm,d )) and beyond ladder (0(G 2 M D )) approximation in the chiral limit. Its 
zeros correspond to the squared bound state masses. Scalar diquarks are seen to be confined 
in O(G 2 0+ ). 

The ansatze of the diquark amplitudes in Dirac space are given by <f>o+(-P) — (P)C 
and <f>i-(P) ~ & P:UJ (P)C, respectively^- When inserting these amplitudes into eq. (|T9D, 
the actual solution of the Bethe-Salpeter equation is obtained in the same manner as for 
the mesons. 


4 Results and Discussions 

In this section we present and discuss our numerical results. The bound state masses for 
the mesons and diquarks in ladder (0(Gm,d)) and beyond ladder (0(G 2 M D )) approxima¬ 
tion are shown in table (|2|). As indicated in the previous section, Goldstone’s theorem is 
manifest for the pion; a massless bound state is found in the pseudoscalar channel when 
the gap equation and the Bethe-Salpeter equation are consistently solved in the same 
order of G = G n (fig. @). When a small current mass is considered in the gap equation, 
i.e. chiral symmetry is explicitly broken, a finite pion mass is obtained. A comparison 
of the pion mass in order G n with the one obtained in order G 2 , entails slightly more 

4 C denotes the charge conjugation matrix, obeying Dy^C -1 = — yj. In the Dirac representation it is 
given by C = *y 2 y°. 


10 








Solution of Vectormeson/1 "-Diquark BSE 



Figure 4: The inverse of the propagators for the vector meson (1 + ) and the axialvector 
diquark ( 1 “) is shown in ladder ( 0(Gm,d )) and beyond ladder [O^Gf^j^]) approximation. 


attraction beyond ladder approximation. 

For scalar diquarks we obtain the conjectured confinement effect: While to order G 0 + 
bound state masses slightly below the pseudothresholds are found, the large repulsive term 
in the Bethe-Salpeter equation (|l^) removes the zeros from the inverse of the propagators 
and therefore confines the scalar diquarks (see fig. ([|)). This behaviour is found in the 
chiral limit as well as in the calculation where a finite current quark mass is considered. 
We observe furthermore, that the inclusion of the subdominant Bethe-Salpeter amplitude 
(~ 75 If C) is necessary to obtain this feature. Working only with the dominant amplitude 
leads in 0{G^+) to scalar diquark masses which, however, lies above the pseudothreshold. 

While the pion (and also the scalar diquark) properties are very much dictated by 
chiral symmetry constraints, the situation is somewhat ambiguous in the vector meson 
and axialvector diquark channel because there is no symmetry available acting as a guiding 
principle for detecting the relevant structures. In this respect, we find that the results 
depend on details of the model, e.g. the ansatz of the amplitude eqs. (|b|,[T^) (one might 
also think of using only as the relevant Dirac structure of the amplitude), the values of 
the coupling constants and even worse, the details of the regularization scheme (see the 
remark in Appendix A). Our choice of G PtU} = G and G i- = 2 G (leading to equal masses 
of vector mesons and axialvector diquarks in ladder approximation) and an ansatz of the 
amplitudes (|T4|), which is ab initio transversal to the total momentum allows at least a 
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Bound state mass in 0(Gm,d ) 

Bound state mass in 0(G 2 MD ) 

Quark Mass 

M = 0.3873 

M = 0.40 

M = 0.4194 

M = 0.4295 

Pion 

0 

0.207 

0 

0.182 

0+ Diquark 

0.762 

0.786 

unbound 

unbound 

p/u; 

0.612 

0.630 

0.737 

0.749 

1~ Diquark 

0.612 

0.630 

unbound 

unbound 


Table 2: In this table we show the bound state masses of the mesons and diquarks in 
0 (Gm,d ) and 0(G 2 M D ). We use the parameters indicated in table ([/]) and values of the 
coupling constants as explained in the text. The infrared cutoff is fixed to p = O.lGeV. 
For the diquarks “unbound” means, that the inverse of the propagators has no zero below 
and slightly above the pseudothreshold P 2 = 4 M 2 , see figure ([|), (0). 

qualitatively comparison of vector meson and axialvector diquark properties in 0(G 2 M D ). 

We then find that the vector meson mass in 0{G 2 pu f) is about 100 MeV larger than 
the mass in 0(G P)U ). For axialvector diquarks the inverse of the propagators has no zero 
below the pseudothreshold (see fig. 0), but the confining effect is not so clearly seen as 
in the scalar diquark channel, where the inverse of the propagator (as a function of P 2 ) 
to 0 (Gq+) will obviously never cross the P 2 — axis, see fig. Here, however, it may be 
possible, that far above the pseudothreshold a zero of the inverse of the propagator signals 
a bound diquark state. On the other hand, as mentioned in sect. (2), the model should 
not be used in the region P 2 4M 2 . Nevertheless, also in the axialvector diquark channel 

the repulsive term in the diquark Bethe-Salpeter equation has a significant effect, which 
may eliminate the axialvector diquark poles. Due to the absence of tensor mesons (and 
therefore also tensor diquarks) in the NJL model there is no a priori subdominant Bethe- 
Salpeter amplitude for axialvector diquarks which could be included in the calculation 
and which may amplify the confining effect in analogy to the scalar diquark case. 


5 Summary and Conclusions 


In this paper we investigated, in the context of a NJL model with an infrared cutoff, the 
conjecture of ref. |T(J that diquarks may be confined in a truncation scheme, which is 
beyond rainbow-ladder approximation. For scalar diquarks we clearly found this feature, 
due to the color structure a large repulsive part of the Bethe-Salpeter kernel eliminates 
diquark bound states, whereas in the pion channel, due to different color factors the same 
term is compensated by attractive parts of the kernel. Vector meson masses are moder¬ 
ately shifted to higher values relatively to their masses obtained in ladder approximation. 
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For axialvector diquarks it is not possible to answer the question of confinement con¬ 
clusively, but the observed behaviour also signals severe changes of the inverse diquark 
propagator in this channel. To conclude, the proposed confinement mechanism for di¬ 
quarks also work in the NJL model in these channels where the model is suitable: the 
pion and scalar diquark channel. 

It will be interesting to see, if this mechanism of cliquark confinement can be confirmed 
in a calculation including a realistic gluon propagator, thereby “interpolating” between 
the NJL and the Munczek-Nemirovsky model. 

Although the question how and whether diquarks are confined is an interesting que¬ 
stion by itself, the outcome of such an investigation has to be included in a baryon 
calculation, where diquarks enter as effective constituents which are preferably confined. 
Together with confined quarks such a description leads to baryons which have no quark- 
diquark thresholds any more |TI§. In such an effective diquark-quark model, besides the 
diquark propagator also the cliquark Bethe-Salpeter amplitude enters, which can also be 
extracted from a diquark calculation presented here or the one of ref. [TO . 


Acknowledgement We thank L. v. Smekal for helpful discussions and critical reading 
of the manuscript. Useful comments by H. Weigel are gratefully acknowledged. 

A Proper-Time Regularization with an Infrared Cut¬ 
off 


In order to explicitly demonstrate how the quark confinement mechanism via an infrared 
cutoff [TTJ works, we show here the solution of the pion Bethe-Salpeter equation in detail. 
For simplicity we will restrict ourselves to the ladder approximation and consider only the 
dominant Bethe-Salpeter amplitude in eq. (0. The Bethe-Salpeter equation ( jT2| ) then 

reads 


t / j—41 . r d 4 k 

>MP)=3 


7 „S(k + iP)7 5 S(i - 


( 20 ) 


Applying the standard techniques for the evaluation of one loop integrals leads, after a 
Wick rotation, to the expression 




/■' r d A k E T 1 (1 - x)P 2 

lo C X J (2^ [ k 2 E + M 2 + (fc| + y 2)2 




Y 2 = x(x — 1)P 2 + M 2 . 


( 21 ) 

( 22 ) 
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Regularization is performed by using the prescription ( 0 ) to obtain 


Via(P) = ^4 


i'dx Pis lpi 

/» 4 J (2 tt)- 


exp(-s(i| + M 2 ))+ 


1 - i)P 2 ) sexp(-s(fc| + V' 2 )) 


( 23 ) 


Note that p 2 = 0 is the usual NJL model expression in proper-time regularization. I 11 eq. 
(p3|) the integration over the loop momentum k and over the proper-time s can be done 
analytically and the final expression (the integration over the Feynman parameter x must 
be done numerically) for the pion Bethe-Salpeter equation is given byQ 




I 67 r 2 3 

p2 

~ 2 AP Jo 


-4 M 2 G* 


M 2 . M 2 n 

r(—1, -j^) - r(—1, -j-p) 


1 / V'2 y 2 

dxim^-no.— 


(p)- 


(24) 


Here the s— integrals provided the expressions involving the incomplete Gamma function, 
e.g., r( 0 , Y 2 / A 2 ) which becomes complex for momenta P 2 larger than the threshold: The 
solutions of Y 2 = 0 for a fixed P 2 are given by 


1 1 M 2 

11/2 ~ 2 ± V 3 - "p 3 "’ 


(25) 


i.e. below the threshold at P 2 = AM 2 they are complex and therefore outside of the 
interval [0,1], the integration region of the Feynman parameter x. Above threshold, 
however, we find Y 2 < 0, for 0 < X\ < x < X 2 < 1, thereby leading to an imaginary part 
of the incomplete Gamma function [ 2H| . 

Y 2 /- 

lmr(0, —) = ±i7r\x2 — Xi\ — ±.mJ ( P 2 — AM 2 ') IP 2 , (26) 

A z v 


displaying the well-known square-root behaviour at the threshold. Since the imaginary 
part is proportional to \x 2 — X\\ but independent of A, it is quite obvious that in the 
expression 



(27) 


which appears in eq. fl24|) the imaginary part arising for P 2 > AM 2 is canceled if /i is kept 
nonzero. 

Although the introduction of an infrared cutoff within the proper-time regularization 
scheme seems to be a convenient way to avoid unphysical quark thresholds in correlation 

5 When comparing this expression with the first line of the gap equation (||) in the chiral limit one 
observes the manifestation of Goldstone’s theorem. 
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functions, one must be very careful when applying this scheme in general. In ref. [^lj it 
has been observed that T(0, Y 2 / A 2 ) has additional unphysical singularities in the complex 
P 2 —plane (which is quite obvious, see [^0|), located on the imaginary P 2 —axis far beyond 
the applicability region of the model, |P 2 | < A 2 . Since their location depends on A 2 , 
they are not canceled by an infrared cutoff. The expression r(0,U 2 //i 2 ) produces even 
more unphysical singularities in the correlation function, also located on the imaginary 
P 2 —axis. Although these singularities play no role if only the Bethe-Salpeter equation 
is solved, they certainly will affect e.g. for m f actor calculations, where a second external 
momentum can shift the singularities into the physically relevant region of the complex 
plane. This indicates, that the proper-time regularization with an infrared cutoff must be 
taken with great care. 

Furthermore, we want to add the following remark: Although the proper-time reg¬ 
ularization scheme is an unambiguous prescription when applied within a bosonization 
approach |[22|| where the fermion determinant is regularized, there are certain ambiguities 
if the method is used on the level of Feynman diagrams 6 . These ambiguities are connected 
with various momentum shifts during the calculation. Since this problem does not arise 
in the gap equation (^), in the pion (and scalar diquark) channel one usually does the 
calculation in a way, which leads to an analytic manifestation of Goldstone’s theorem. 
However, in the vector meson (and axialvector diquark) channel, chiral symmetry can 
not be used to resolve the ambiguities. Therefore the results in these channels are not 
unambiguous. 


B counting beyond ladder approximation 


In the meson and diquark Bethe-Salpeter equations different color factors arise which in 
turn lead to the differences in meson and diquark channels. In color space, the meson 
Bethe-Salpeter flT2|) equation can be denoted by 


<P?k = + ifPi )^LCj*LC + (Vs )K^mCn^n^fci) ( 28 ) 


where the (t a ) are the generators of color SU(IV C ). The meson amplitude <p^,(A = 
L, 1, 2, 3, where L denotes the ladder part and 1, 2, 3 the three terms beyond ladder ap¬ 
proximation), describing a color singlet state, is given by the ansatz tpfj = (p M 5ij / \fN~ c - 
Using this ansatz leads to the following color factors: 

(I) : = \^WT' ° (Ay (29) 

6 The usual bosonization approach is deeply connected with the rainbow-ladder approximation and 
cannot be used in a truncation scheme like the one considered in this paper. 
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( 1 ) 

( 2 ) 

(3) 


1 

iVc 

1 






x- 


4 JV? 


TV 2 


4iV c 2 






N? 


w 


0 ( 1 ) 

0 ( 1 ) 

0 ( 1 ), 


(30) 

(31) 

(32) 


i.e. for N c = 3 the factors appearing in eq. fll2|) are obtained. In the spirit of an 1/N C — ex¬ 
pansion we observe that the term in ladder approximation would be of leading order 
(0(iV c )), while the terms beyond ladder approximation are subleading, (0(1)). 

The diquark Bethe-Salpeter equation (|1^) in color space can be written as 


<P?k = M)lit a iflk + {(Vl'jliti/rnktnmtfn + C/nfcC + (<Py )> ( 33 ) 


where the prescription (^) has been used to derive it formally from the meson Bethe- 
Salpeter equation (note, (t“ ) T = t,X). The projection on color antitriplet diquark states is 
most conveniently done by using the ansatz <p? = y> D e“-/ \/N c — 1 of the diquark ampli¬ 
tude, involving the normalized Levi-Cevita tensor eXj\/N c — 1. Then the resulting color 
factors in the diquark channel are given by 


(L') 

(10 

(20 

(30 



c j.a j.a a 


N c + 1 
2 N c 


zed 



c j.a j.b j.d ±b d _ 

^li^ij^mk^nm^ln^jk 


c ib ±d ib id d _ 

^li^im^mj^nk^ln^jk 


c ib id ib id d _ 

^li^im^mn^nj^lk^jk 


N C + 1 red 

(2 N c y ’ 


N c + 1 
(2 N c y 


(1 + N C 


N c + 1 „ cd 
(2 N c y ’ 


0 ( 1 ) 
0(1/N C ) 
N%) S cd , 0(N C ) 
0(1/N C ). 


(34) 

(35) 

(36) 

(37) 


For the case N c = 3 we obtain the color factors appearing in eq.(|T9l). As mentioned in 
section(^D, the different terms in the diquark channel shows a rather interesting N c —► oo 
behaviour. In leading order one would just keep the term (2') which is of 0(N C ). However, 
the whole contribution of this term in the diquark Bethe-Salpeter equation is repulsive 
and thereby eliminates diquarks totally. The attractive ladder part of the Bethe-Salpeter 
equation contributes in this scheme only in subleading order, 0(1), and the attractive 
terms beyond ladder approximation only in order 0(1/N C ). 
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